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Abstract 

The study of quenched random systems is facilitated by the idea that the ensemble averages 
describe the thermal averages for any specific realization of the couplings, provided the system is 
large enough. Careful examination suggests that this idea might have a follow, when the correlation 
length becomes of the order of the size of the system. We find certain bound quantities are not 
self-averaging when the correlation length becomes of the order of the size of the system. This 
suggests that the strength of self-averaging, expressed in terms of properly chosen signal to noise 
ratios, may serve to identify phase boundaries. This is demonstrated by using such signal to noise 
ratios to identify the boundary of the ferromagnetic phase and compare the findings with more 
traditional measures. 

PACS numbers: 05.50.+q, 64.60.Cn, 75.10.Nr, 75.10.Hk 



Quenched random systems like random field systems, random bond systems, spin glasses 
etc. are known to be extremely difficult, because of the necessity to perform quenched 
averages. In fact, it has to be realized, of course, that this necessity, is not of physical origin 
but rather of mathematical convenience, since a given chunk of matter that is measured 
experimentally has a given single realization of the disorder. The idea of self-averaging 
allowing ensemble averaging is based on the fact that we are dealing with very large systems 
It is argued that the system can be broken up into subsystems large enough to be 
considered independent of each other. This basic assumption indicates that phenomena 
connected with breakdown of self-averaging will appear, as the correlation length will become 
of the order of the linear size of the system. Indeed, such behavior was first observed by 
Dayan et al Q], using a technique first suggested by Berker and Ostlund jjj. The subject 



of breakdown of sel 



phenomenon 
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-averaging was treated later in a number of papers devoted just to that 



Obviously, lack of self- averaging results in severe difficulties but 
the fact that it is connected with the divergence of the correlation length suggests that it 
can provide, perhaps, an independent measure to distinguish between the disordered and 
the ordered phase. In this article we show that this is indeed the case. 

To be specific consider the random field Ising model described by the Hamiltonian 

H = -J) j ViO-j - KiJi. (1) 

<i,j> i 

The pair < i,j > denotes a nearest-neighbor pair on a cubic lattice. The field configuration 
is assumed to be governed by the distribution, 



PW =n^-p(-^ 



(2) 



We take J and h to have the dimensions of energy. 

Consider now some quantity a and its ensemble average, ag. Define next the variance of 
a within the ensemble, 

a a = {[a 2 ] -a E 2 }^ (3) 

where [. . .] denotes ensemble average. The parameter 7 describing the strength of self- 
averaging in the system (or rather its weakness) is defined as the signal to noise ratio 

la = • (4) 

(Try 
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It is clear that for the disordered phase 7 is infinite in the infinite system. Our former dis- 
cussion suggests that possible breakdown of self-averaging may be traced back to situations 
where the correlation length, £, becomes of the order of the linear size of the system. The 
usual state of affairs is that the correlation length diverges just on the boundary between 
the ordered and disordered phase and consequently any break down of self-averaging may be 
observed only in the vicinity of the boundary 4], ^ |(| . In the random field case the situation 
is quite different. The correlation length is of the linear size of the system everywhere in the 
ordered phase and not just at the boundary of the phase Q]. 

To understand the possible different behaviors of 7 consider first the case where a is 
the free energy. Suppose we break the system into mutually non-interacting sub systems of 
linear size L'<L and evaluate the free energy of the system as the sum of the free energies 
of the sub systems. This generates a relative error in the free energy that is of the order of 
1/L' regardless of the size of the correlation length. This suggests that nothing spectacular 
will happen to the signal to noise ratio of the free energy when we cross from the ordered 
to the disordered phase. To see if this is true we need to know, first, where is the ordered 
phase boundary. In Fig. ^1 we give the logarithm of the susceptibility. The ordered phase 
is presented as an upper plateau. In Fig. we give the Laplacian of the logarithm of the 
susceptibility with respect to the coordinates given in Fig. ^1. The white line drown in Fig. 
\I]p is where the absolute value of the Laplacian, which is related to the curvature, is maximal 
and is viewed as the phase boundary. Our result for the boundary is consistent with phase 
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diagrams (or usually parts of it) obtained by others 

In Fig. El we give the signal to noise ratio of the free energy. Nothing spectacular happens 
in the vicinity of the phase boundary. The signal to noise ratio in the susceptibility is another 
matter. It remains small anywhere in the ordered phase. For the sake of completeness we 
give in Fig. |3] the temperature dependence of the signal to noise ratio associated with 
the susceptibility as obtained from the results of Ref. j2|. If only signal to noise ratios in 
divergent quantities showed some interesting behavior it would not have been very interesting 
as indeed it would have provided no additional information to that obtained by calculating 
the average itself. It is thus of interest to consider physical quantities that do not diverge 
in the ordered phase or on its boundary. Our choice is motivated by our previous work on 
real space renormalization of the three dimensional random field Ising model (RFIM) [18]. 
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FIG. 1: The logarithm of the susceptibility (a), and its Laplacian (b) everywhere, in dimensionless units, 
as a function of the dimensionless temperature, ksT/ J, and the dimensionless strength of the random field, 
h/J. In (b), a pseudo-critical line of the finite 32 x 32 x 32 system is represented by the white line marking 
the dark ridge separating the two brighter plateaus. 
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FIG. 2: The signal to noise ratio of the free energy of a system of size 32 x 32 x 32. Everywhere it is of 
the order of 10 4 . 

In that article we have obtained, for the first time, the "connected" spin-spin correlation, 

T E (r) ee [r(r)] = [(a(O)a(r)) - (a(0))(a(r))]. (5) 

In Fig. EJwe present Te for the case where the strength of the dimensionless random field 
is h/J = 1, for linear sizes of L = 64 and L = 128 and three different temperatures. It is 
clear that as temperature is lowered the noisiness increases up to the lowest temperature 
(Fig. 0k), where even the trend of decay as a function of r cannot be observed any more. 
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FIG. 3: The signal to noise ratio in the susceptibility for a system of 32 x 32 x 32 and for a given 
ratio of ksT/h = 2, obtained from data presented in Ref. [2j. The critical dimcnsionless temperature is 
k B T/J = 4.2. 

Since Te must be a decaying function of the distance this implies that self-averaging must 
be broken in Fig. 

An attempt to measure T(r) can be made, perhaps, along the following lines. An ex- 
periment is carried out on a given chunk of matter with some fixed position dependent in- 
teractions, so that the measured quantities are spatial averages on that system rather than 



ensemble averages. The spin-spin q dependent corre 



ation function, S(q), can be obtained 



experimentally from neutron scattering o, o, q, y, q> making the proper allowance 

that the measurement is done on a dilute anti-ferromagnet in an external constant field 
rather than on a ferromagnet in the presence of a random magnetic field. The question 
now is whether the "connected" correlation, yftf h can be unpicked from S(q). To do that a 



relation between the two has to be assumed 



22|, |2^J. The spatial rather than the ensemble 



average of the Fs(r) is given by the Fourier transform of x(<?)) 

rs(r) = J dqx(q) exp(zg • r). (6) 

The evaluation of r,s(r) has to be repeated for a large enough number of systems in order 
to obtain reasonable ensemble averages. We will concentrate in the following on the sim- 
pler quantity r(0) and make our point by obtaining the corresponding 7 from real space 
renormalization of the random field Hamiltonian ((TJ. First, however, let us obtain a crude 
picture of what to expect. Consider some local bounded quantity, Tj, such that the ensemble 
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FIG. 4: (Taken from Ref. e lSj) The average spin-spin correlation function, Te, is shown as a function 
of the distance, r, taken along the main axes of the lattice and measured in units of lattice constant. The 
external field, H, is zero while the standard deviation of the RF is h/J = 1. As indicated on each figure, 
the two different levels of gray scaling correspond to systems of different linear size, L, indicating the size 
independence of Te for these temperatures. The points of the larger system, with L = 128, appear to be 
more scattered since it is averaged only over 3500 realizations, while the smaller, L = 64, system is averaged 
over 10000 realizations. The different figures correspond to different temperatures. Note the broadening of 
as the temperature is reduced [(a) and (b)] towards entering the ordered phase, at about ksT/J = 3.8, 
untill, due to noise [as in (c)], no trend can be detected. 



average of Tj is T^O). Consider next the spatial average 



(7) 



where N is the total number of sites in the system. As long as self-averaging is not broken 
Ts{0) is expected to "equal" its ensemble average, T E (0). (By "equal" we actually mean 
that the corresponding signal to noise will be very large.) We can get an idea of the variance 



in Ts(0) by noting that the sum in the numerator on the right hand side of Eq. (JJJ) is a sum 
of correlated random numbers that are correlated over a distance £. This implies that the 
corresponding signal to noise, 7 r , is given by 

7r ~ (£) 2 (8) 



for a d dimensional system. (Indeed, when approaching the boundary between phases from 
the disordered phase a more accurate exponent should be used reflecting the nature of 
algebraic decay of correlations within the correlation length. This will replace d/2 in Eq. 
(JHJ) by another (smaller) positive power of £.) In any case the signal to noise must decrease 
considerably in the ordered phase compared to its values in the disordered phase. (Note that 
the free energy, cannot be written in the form given by Eq. (J2J) and therefore the argument 
given above does not apply to signal to noise ratio of the free energy.) It is clear ,however, 
from the argument above that relation (JHJ) is of a generic nature and does not depend on the 
specific quantity used. The ratio of two signal to noise ratios or two variances corresponding 
to two different quantities will thus not be essentially different in the two phases. 

We describe next our procedure for obtaining the signal to noise ratio of T(0) numerically, 
which is based on Casher-Schwartz RSRG j^j]. As other real space techniques, it provides 
simple, one step, recursion relations for translational invariant systems that enable the ex- 
traction of the critical exponents. For random systems, the recursion relations obtained, 
using any renormalization scheme, involve the distribution of couplings or, equivalently, all 
the parameters defining it (e.g. moments, correlations etc.). In this approach, the recur- 
sion relations are truncated to obtain relations involving only the mean and the variance and 
keeping the random couplings independent [^Ej 21, 2s| - The method suggested by Berker 
and Ostlund [3] overcomes the most difficult problem arising in the approach described above 
of ignoring the correlations generated by the renormalization procedure (or projecting them 
on the variance A given realization is chosen on a finite system. Renormalization is 

then used to reduce the size of the system to a size where brute force calculation is possible. 
(The fact that a lot can be learned from considering specific realizations is also stressed in 
a recent paper on the RFIM by Wu and Machta [29].) After the required thermal averages 
are obtained for a given configuration of the disorder, the ensemble average is obtained by 
repeating the procedure for many configurations and averaging. Only thermal averages of 
functions of the spins surviving the renormalization seem to be obtainable directly. Berker 



and coworkers 



3 it |32j were using, however, the chain rule to approximately recover 
thermodynamic densities of the original system from the renormalized couplings of the re- 
duced system. This enables to obtain averages which involve not only surviving spins, but is 
limited to thermal averages of spin products which appear in the original Hamiltonian (such 
as nearest-neighbor pair products). In the following we will see how averages involving non 
surviving spins such as T(r) can be calculated [l^ |. 

We start with a set of N = L 3 Ising spins, with L = 2 n , situated on a three dimensional 
cubic lattice with periodic boundary conditions. We generate a realization of the distribution 
(J2J) and then perform the Casher-Schwartz procedure (2^ n — 1 times. At each step of the 
renormalization, the lattice remains cubic while its linear size is reduced by a factor of two. 
The couplings remain nearest-neighbor but become position dependent and each spin is 
affected by a renormalized position dependent field. Finally, we are left with a system of 
2x2x2 spins on which we perform direct calculations. Consequently we obtain 8 (c^'s, 
corresponding to the spins that survive the renormalization procedure. We define a local, 
position dependent, quantity 1^(0), 

r *(°) = W {(Tl)K (9) 

where hi is the original field on the site i. The actual ensemble averages are performed by 
repeating the calculations for many realizations of the randomness, summing up the results 
and dividing by the number of realizations. The statistics is improved by using, instead of 
the local quantity defined above, its average over the spins surviving the renormalization 
procedure, 

r (°) = ^E^> ( 10 ) 

where i runs over the 8 surviving spins. (Note that the ensemble average of T(0) is indeed 



the "connected" self correlation r^(0) \Yt 



The variance is 



^r(o) = {[r 2 (o)]-r s 2 (o)}^. (11) 

The signal to noise in T(0) is presented in Fig. E3 The strong reduction in the ordered 
phase is very clear. We proceed now to check relation (jHJ), which is supposed to be generic. 
We define first the quantity 

^86 

r (i) = 4^£ EMM. (12) 

H i=l n(i)=l 
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FIG. 5: The signal to noise ratio in 7r(o) ( a ) an d its Laplacian (b) everywhere. The strong reduction in 
the ordered phase, shown in (a), is evident. As in Fig. ^ in (b), the white line marks the maximum of the 
Laplacian and taken as the phase boundary. 




FIG. 6: The ratio of ar(o) an< ^ °r(i)- This does not show any pronounced dependence on the correlation 
length. 

where n(i) runs over the six sites that are nearest-neighbors of i in the original lattice (all 
sites that, in the original lattice, are at a distance r = 1 from i). In Fig. El we depict the 
ratio between 0T(o), the variance of r(0), and 0r(i), the variance of T(l). This ratio does not 
show any spectacular behavior when the transition line is crossed. 

The conclusion is that although quenched random systems are not self-averaging every- 
where j^l 0, S, 0, , this undesirable feature can be turned into a useful tool. The strength of 
self- averaging expressed in terms of signal to noise ratios in quantities that in pure systems 
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may not even be interesting, becomes a useful quantity to study. These supply additional 
measures of the order in the system, because, local bounded quantities lead to signal to 
noise ratios that depend on the correlation length in a generic way. 
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